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Abstract 

Two numbers m and n are considered amicable if the sum of their proper divisors, 
s(n) and s(m), satisfy s(n) = m and s(m) = n. In 1981, Pomerance showed that 
the sum of the reciprocals of all such numbers, P, is a constant. We obtain both a 
lower and an upper bound on the value of P. 



1. Introduction 

Since at least the time of the Ancient Greeks, amicable numbers have enjoyed the 
attention of mathematicians. Let s(n) denote the sum of the proper divisors of 
n, that is, s(n) = o~(n) — n, where <j(n) denotes Euler's sum-of-divisors function. 
Then an amicable pair is a pair of distinct integers m, n, such that s(n) = m and 
s(m) = n. We will also refer to any integer which is a member of an amicable 
pair as an amicable number. The smallest amicable pair, (220,284), was known to 
Pythagoras c. 500 BCE. The study of amicable pairs was a topic arising often in 
Medieval Islam; as early as the 9th Century, Thabit had discovered three pairs, 
including (17296, 18416) - a pair which was rediscovered independently by Borho, 
Ibn al-Banna', Kamaladdin FarisI, and Pierre de Fermat [I]. In the 18th century, 
Euler [T] famously advanced the theory of amicable numbers by giving a table of 30 
new pairs in 1747 (one of which is, in fact, an error - see P2]), and a larger table 
of 61 pairs, together with a method of generating them, in 1750 [2]. 

Today, although much is known about amicable numbers (and their less popular 
friends, sociable numbers), there is still a lot that we don't know. For example, 
although the best upper bound on their density was given by Pomerance |llj in 
1981 (see below), no known lower bound exists - indeed, it has not been proven 
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that there are infinitely many. Considerable work has been done on questions such 
as the properties of even-even pairs, odd-odd pairs, pairs for which each number 
has exactly two prime factors not contained in its pair, and many more complex 
variations. For a nice survey on amicable numbers, see [3]. One interesting fact, 
which motivated the present work, is that the sum of the reciprocals of the amicable 
numbers converges, that is, 



^r* a 

a amicable 



= P<oo. (1) 



This is a consequence of the bound shown 28 years ago by Pomerance [TTJ , but to 
date no bounds have been given on the value of this sum. In this paper, we provide 
an upper bound and a lower bound for the value of this sum. 

1.1. The distribution of amicable numbers 

Following [TT], we define A(x) to be the count of amicable numbers not greater 
than x. In 1955, Erdos [3J showed that the amicable numbers have density zero. 
Pomerance showed [TU] 

A(x) < a; exp ^— cy/log log log x log log log log x 

for a positive constant c and all sufficiently large x. Then, in 1981, Pomerance [UJ 
improved this result to 

A(x) < C X (2) 

exp ( (log a;) J 

for some constant C and all sufficiently large x. From this and partial summation, 
it is clear that we have ([1]). Since Pomerance was the first to show this sum to 
converge, we refer to its value, P, as Pomerance 's constant. 

Having established that P is finite, it is natural to ask about its value. We prove 
the following theorem. 

Theorem 1.1. Pomerance's Constant, P, the sum of the the reciprocals of the 
amicable numbers, can be bounded as 

.0119841556 < P < 6.56 x 10 8 . 

It should be noted that the upper bound in this result can be improved by a 
more careful choice of the functions in Table [JJ on page [3] It is unclear how to fully 
optimize this argument, but the authors have been able to show P < 3.4 x 10 6 . 

To establish Theorem ll.il we split the sum defining P into three parts: 



E 



E 1 + E 1+ E -• ( 3 ) 

n * — ' n £ — ' n * — ' n 

amicable n amicable n amicable n amicable n 

n<10 14 10 14 <n<cxp(10 6 ) n>exp(10 6 ) 
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We evaluate the first sum directly in section [2] to establish the lower bound for P. 
In section [31 we bound the second sum. We then modify the argument in [TT] to 
address the final sum above in section [5] 

Throughout, p, q, and r will denote primes, P(n) the largest prime factor of n, 
and ip(x, y) the number of y-smooth numbers up to x, that is, the size of the set 

{n < x\P(n) < y}. 

Finally, a and b will represent amicable numbers. 

We will repeatedly make use of the following functions. These will be redefined 
at the appropriate place, but the reader can use this table as a helpful reference. 

x = exp(10 6 ) 

y = exp(26000) 

c = c(x) = 1 — (logx) -1 / 6 (log log x) -1 

c = c(10 6 ) > 0.99276 

g =a{y) = 1-1/(2 logy) 

= l(x) = exp ^(logx) 1 ^ 6 ^ 
= L(x) = exp (0.1882 (log x) 2/3 log log x 



i 



Table 1: A table of functions used in the proof of Theorem ll.il 



2. A Lower Bound Some Reciprocal Sums 

Determining a lower bound on P is straightforward - we need only sum the recip- 
rocals of any subset of the set of amicable numbers to find one. Let Pj be the sum 
of all integers not greater than 10 J which are members of an amicable pair, i.e., 

a<10 J 
a amicable 

The current record for exhaustive searches for amicable numbers is 10 14 [9]. In 
Tabled we give values for Pj for various powers of 10 up to this bound. 
From this table, we immediately have 

P > Pu > .0119841556796931142. (4) 

We note that this series seems to converge rather quickly. Indeed, for any j < 13 
with Pj > 0, the difference between between Pu and Pj is less than the difference 



INTEGERS: 10 (2010) 



4 
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Pj 


1 





2 





3 


0.0080665813060179257 


4 


A H1 1 1 A A^l A*7 A A Cid 

O.Ulllb 1 1 zb 14424 (44bO 


5 


0.0117423756996823562 


6 


0.0119304720866743157 


7 


0.0119714208511438135 


8 


0.0119812212551025145 


9 


0.0119834313702743716 


10 


0.0119839922963130553 


11 


0.0119841199294457703 


12 


0.0119841486963721084 


13 


0.0119841542458770555 


14 


0.0119841556796931142 



Table 2: Sums of reciprocals of amicable numbers to 10- 7 . 

between Pj and Pj-i- With this observation, we conjecture that the true value of 
P can be estimated by 

\P - Pu| < Pu - Pis ~ .0000000006338, 

and thus: 

Conjecture 2.1. Pomerance's constant, P, satisfies 

P < .0119841563134. 

This, however, is merely conjecture. We turn now to establish the upper bound 
of Theorem O 



3. The Middle Range 

The amicable numbers in the range [10 14 , exp(10 6 )] are too large to be found ex- 
plicitly, and too small to be amenable to the results we make use of later in this 
paper. While there are some ways to restrict the sum over these numbers, their 
final contribution to Pomerance's constant is small, and we here simply make use 
of the trivial bound: 

y -< io 6 . 

10 14 <n<exp(10 6 ) 
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4. Preliminaries 

We will use a number of explicit formulas of prime functions. For instance, we will 
use the fact (see [HI Theorem 5]) that 

V - < log log x + B H - ? 

for x > 286, where £? = .26149721 . . .. Recall that we have chosen xq = exp(10 6 ). 
This gives that, for x > xq, 

V- < log log a; +.2615. (5) 

p<:e 

We will use a bound on the sum of reciprocals of primes in a particular residue 
class modulo p, where p > 10 . 

Lemma 4.1. For any value of y > 10 14 and any prime p > 10 14 , the following 
holds uniformly: 

1 < 4 + 31og logy 

^ q ~ p 

q= — l mod p 

Proof. We begin with the Brun-Titchmarsh inequality of [8]; namely, for coprime 
integers k and n, the number of primes q < y with q = fcmod n, denoted 7r(y; n, k), 
satisfies 

2y 

n(y; n, k) < — — — — - for y > n. 
tp[n) log(y/n) 

For prime p > 10 14 and (k,p) = 1, we can use the fact that ip(p) = p — 1 to see that 

f M <- 2 y / 2.0001y 

^(y; P= fc ) < TTi — T~T\ — ~~ i — TTT' ( 6 ) 
WJ log(y/p) p\og(y/p) 

Note that the first prime g = —1 mod p is at least 2p — 1, so we need only consider 
p < 2±1. Thus, 2 > 2 - i > 1.999999 for p > 10 14 . We use partial summation to 
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obtain 



^ q ' V J 2p -i t 2 



q<y 

q= — l mod 



2.0001 [ v 2.0001 
< -: ^-r-r + / : t-tt dt 



plog(y/p) J 2 p-lPtlog(t/p) 
2.0001 2.0001 f y dt 



p\og(y/p) p J 2 p-itlog(t/p) 
2.0001 2.0001 , , , , . 

n — r~r\ + log log t/w 

p\og{y/p) p 2p _ 1 

2-0001 / 1 , , , . , , , / 1 

+ loglog(j//p) - log log 2 - - 



p V lo g(z//p) 

Since iogfa/p) ~ lo S lo g ( 2 ~ p) < i- 8093 and loglog(y/p) < log logy, we see that 



E 



1 4 + 31oglogy 



- < 



q p 

q<y 

q= — l mod p 

which proves the lemma. □ 

We will need the following lemma to bound the number of amicable numbers 
up to x which are also y-smooth for some y > yo = exp (26000). Recall that 
a = 1-1/(2 logy). 

Lemma 4.2. For y > yo and a as above, we have 

Iifi-i) <7.65151ogy. 

p<y v y 7 

Before proving Lemma |4~21 we establish a few necessary lemmas. 
Lemma 4.3. For y > yo, 

logy - 1.3325. 

t?i P 

p<y 

Proof. We recall a result from Rosser and Schoenfeld (see [12l 3.23]). Let E — 
-7 - E p (l°gp)/P n < -1-332582275. For y > 319, 

Elogp , _, 1 
<logy + £+- . 
P <y p 2l0 ^y 



Setting y = y , Lemma 1431 follows. □ 
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Lemma 4.4. Let y > yg. Then, for a as above, 

V — < log log y + 1.0859. 

p<y 

Proof. We first note the helpful fact that 

f =pcxp(-logp/(2 logy)) >pexp (-1/2). (7) 

The idea in this lemma is to bound the difference between 1/p and l/p a ■ Let 
f(y) = P V j then we want to estimate f(a) — /(l). Since f'(y) = — (logp)/(y), it 
follows from the mean value theorem that f(a) — /(l) = (1 — a) (log p)/p v for some 
y between a and 1. Since a = 1 — 1/(2 logy), this gives 

/(.)-/(!)< 1 



2 logy V P 7 
e 1 / 2 /logp 



2 logy V P 

by ®. 

Hence, by Lemma 14.31 we have 



( 1 1\ eV2 v^logP e 1 / 2 Aogy- 1.3325\ 

V < — > — — < — — — < 0.8244. 

f^ y \P a P/"21ogy^ p 2 V logl/ J 



So this means that 



V-<l + T- < log logy + 1.0859, 



p u ^ p 

p<y p<y 

by ©. □ 
Lemma 4.5. For y > yo awd tr as above, we have 

y — - — - < .7734. (8) 

^ p°(p° - 1) 

Proof. Since y > yo, it follows that a > .9999. We will bound ((8]) by explicit 
computation and analytic methods. Using 500000 as a useful place to split our 
computation, we write 

V - < V - + V - 

p r ^ ' p<500000 1 yr ' 500000<p<oo y ^ ' 

f°° 1 

< .7733545 + / -— dt 

J 500000 t a \t a — 1) 

< .7733545 + .0000001 

< .7734, 

by explicit computation. □ 
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We may now prove Lemma 14.21 
Proof. 



exp 



exp 



1 



per _ 1) 



< exp (log log y + 1.0859) exp(.7734) 

< 6.4193 logy, 

by Lemmas [43] and 14.41 □ 



We now show that our choice of x > xq is "sufficiently large" to adapt Pomer- 
ance's argument [TT]. Let 

c = 1 — (loga;) -1 ^ 6 (logloga;) -1 . 

For x > xq, it is clear that 

c > c = 1 - (10 6 )~ 1/6 (loglO 6 )" 1 > .99276. 

Lemma 4.6. For any c > c above, we have 



y 1 

^ k c (k c - 1 



, < 1.0225. 

.. (k c - 1) ~ 

fc>2 V ; 

Proof. First, we split the sum into two parts, i.e., 

^ k c (k c - 1) ~~ ^ & c (fc c - 1) + ^ fc c (fc c - 1)' 

fe>2 V ; 2<fc<10 6 V ; fe>10 6 V ; 

For any value of k > 10 6 , we have k c < 1.000002(/c c — 1), so the second sum can be 
bounded by 

1 1 f°° dt 

V — ; < 1.000002 V — < 1.000002 / — < .0000013. 

^ k c (k c -l) ~ ^ k 2c ~ i ln6 , t 2c ~ 

k>w & y ' fc>106 JU> - 1 

A quick computation shows that 

V —— < 1.02247315, 

^ k c (k c - 1) ~ 

and adding these two together completes the proof of the lemma. □ 
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Lemma 4.7. For any c > cq, we have 



y 1 

^ c (p c_ 1) 



< .7877. 

n c ( n c — 1 ) 

p>2 



Proof. As before, we split the sum into two parts, i.e., 

p c f p c _ i) ~~ X! pC (p C _ J) + X] „c foe _ ■ 

Once again, we use an explicit computation to show that 

V — -r- < .7876817684, 

p c fpC _ 1) ~ 

and the second sum can be bounded by 



y 1 



. < .0000013 

„ p c iv c — 1) 
p>10 6 ^ v ^ ' 

as in the previous lemma. Adding these two together completes the proof of the 
lemma. □ 

In section[5j we will use the function £ = exp ^(logx) 1 ^ 6 ^ , which is referenced in 
the following lemma. 

Lemma 4.8. For any x > x and c = 1 — (log log a;) 1 (log a;) X ^ , we have 

V — c < i log log a; + 2.0346. 

p<£ 4 r 

Proof. First, note that 

pC = p i-(io S \o gx )-Hio S xr 1/6 =p exp ( l ° gP ^ 



> pexp 

= p exp 

> .7486p 



(log a;) 1/6 (log log x) j 

4 (log x) 1/6 
(log x) 1/6 (log log a;), 
4 



(log log a;) 



since p < £ 4 and x > e 10 . Then, note that the mean value theorem says that, with 
f(y)=p- y and f(y) = -(Iogp)/p*, 

/(c)-/(l)<^#^- 
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Applying both of these facts, we see that 



y: (- c --)< i.3358(i - C )x: 

By Lemma |4~31 this can be bounded by 



(4(logx) 1 

(1 - c)l. 3358 (log^ 4 - 1.3325^ - 



1.3358 (4 (log xY /6 - 1.3325! 



(log x) 1/6 (log log x) 
5.3432 

< 



log log X 
< .3868. 



Then, we have 



1 1 



t. — < pC Z — < p Z — < \pc p 

p<£ 4 p<^ 4 p<£ 4 V 

< log log £ 4 + .3868+ .2615 

< I log log x + 2.0346. 
6 

by ©. □ 

In [TT| , Pomerance uses a function N(s) which, for 

il>(m)= J] (q + l)q a -\ 

q a \\m 

is the number of m with ip(m) < s and P(jp(m)) bounded by a somewhat small 
function of x. His bound applies only to "large x" , so we must verify that our choice 
of x > xq is large enough. Pomerance showed that 

N(s)<2s c Yl (l-&~ c ) _1 . 

P(fe)<C 4 
k>\ 

We now follow making the bounds explicit as we proceed. 

n (1-^)"= n < 9 > 

P(fe)<£ 4 P(fc)<£ 4 V ' 7 

fc>l fc>l 

Since 1 + < < e 4 , we have (for i = pr + + • • • ) 

11 /ll 

l + F + ^ + -"<exp - + -^ 
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Putting this into (J9j) gives 



n (i-^r 1 ^ 



exp 



P(k)<£ 4 
k>l 



E k ~ c + E r 2c + fc ~ 



3c 



. p(k)<e 4 p(k)<e 4 

\ k>l k>l 



(10) 



We may bound 



E (^ 2c + fc ~ 



3c 



P(fc)<£ 4 
fe>l 



• ) < y Trr — tt < 1-0225, 
' ~ ^ k c (k c - 1) ~ 



by Lemma 14.61 Setting aside this secondary term, the remaining term in the expo- 
nential in (fTOj) satisfies 

E k-^nii-p-r 1 - 

P(k)<l 4 p<l 4 
k>l 

Bounding this product by the same argument as in ((Q)) and (fT0|) . we have 

P <e 4 \p<£ 4 p<e 4 ^ ^ ' J 

< exp (.7877) exp ^ p" 

\p<i 4 



by Lemma 14.71 An application of Lemma 14.81 shows that this product is bounded 
by 



JJ (l-p _c ) 1 <expQloglo.. ; .r-2.s22:-!) .■.xp(2.N22:-!) (lo^.r) 



1/6 



Thus, 



J| (l - *T C ) 1 < exp (l.0225 + exp(2.8223) (logs) 



i/(> 



p(fe)<£ 4 
fe>i 

We will need the bound in (j9]) to hold for 

s > L = exp (0.1882 (logx) 2/3 log log: 
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If s > L, the argument above shows that the inequality 

< 2s - 1 exp ( 1.0225 + cxp(2.8223) (logx) 1/6 



= 2 exp ((c - 1) log s + 1.0225 + exp(2.8223) (log x) 1/6 
< 2 exp (-.1882 (loga;) 1/2 + 1.0225 + exp(2.8223) (logx) 1/6 ) 

2e 1.0225 

holds for x such that 

2 + exp(2.8223) < .1882 (logx) 1/3 , (11) 

which, in turn, holds true for x > xq. 

Lastly, we will make use of the fact that for k > 2, 



n 

n>k 



5. Large Amicable Numbers 
5.1. Outline of the Proof 

Pomerance's proof of ([2]) rests on an argument that five different types of numbers 
do not contribute much to the sum defining P, with a careful count of amicable 
numbers among the remaining integers up to x. For the sake of brevity, we include 
only the statement of each assumption, and the resulting bound which must be 
made explicit to bound the value of P. We make a similar argument here, though 
we bound A(x) by a much smaller function so that our bound will apply for smaller 
x. Following the notation of [11] , we will use our previously defined functions 

£ = exp ((logic) 176 ) and L = exp (.1882 (log x) 2/3 log log . (13) 

Note that I is smaller and L is slightly larger than their equivalents in [TT]. It will 
also be useful to have z = 2x log log x. 

Call an integer n admissible if it satisfies each of the five conditions given in |11) . 
We summarize them here. 

(i) The largest prime factor of n and s(n) are both at least L 2 . 

(ii) If k s divides n or s(n) with s > 2, then k s < £ 3 . 

(iii) If p | gcd (n, c(n)), then p < £ 4 . 
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(iv) The integer n satisfies 



n , s(n) 

> L and — -— — > L. 



P(n) ~ P{s(n)) 
(v) If m = -pj~) and m! = p(f "2)) ' then we have P (a(m)) > £ 4 and P (cr(m')) > 

In showing (|2J), Pomerance proved the following (see [11] for details). 
Theorem 5.1 (Pomerance). The set of integers n < x which are not admissi- 
ble is O (x/ exp ^(log xf^^j^ . Specifically, the number of integers n < x failing 
conditions (i)-(v) above are bounded by the functions 
(i) ^(z,L 2 ), 



(a) z y) k 



(m) Y] Y] — , 

, PQ 

p>£ 4 q=~l(p) 
qKx/p 

4xloglogx 

(m) , and 

L 

4e 1 0225 a;loglogx (1 + 2 log logs) ^ 1 
(V) P \ p' 

respectively. 

Unfortunately, in |llj . the bounds in (v) are only valid for very large x. By 
choosing the values of I and L in (fl"3"|) . which are much smaller and slightly larger, 
respectively, than their counterparts in |llj . we can make the argument apply for 
all x > xq — exp (10 6 ). In making these bounds explicit, we bound the quantity in 
(i) by the function C\xj t, the function in (ii) by C2X/I, and so on. 

Let S denote the set of all admissible integers. Pomerance demonstrated that 
the count of amicable numbers in S can be bounded above by 

£ £ £ £ 5«<^ ,») 

r>£ 4 q=-l(r) m=0(q) q' = -l(r) 
q<x m<x q' '<z 

for some constant Cq (here (n) should be taken to mean mod n). 

Then, it is clear that we may take C = C\ + C2 + ■ ■ ■ + C§ in the bound 

A{x) < C j-^- — y (15) 

exp ^(logir) 1/6 J 

for x > xq, and so it remains to find explicit values for each of these constants. 
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5.2. Evaluating the constants 

We address each of the constants given in subsection 15.11 

C\\ We use an explicit version of the method of Rankin to bound ip(z,L 2 ). Given 
that x > Xq, we also have that L 2 > y . Rankin's method is based on the 
observation that, for any <j > 0, x > 1, and y > 2, 



?i>l n>l 
P(n)<q P(n)<q 



X 

p<y 



(16) 



From Lemma T4. 2 1 is follows that 

1>(x, y) < 6.4193 (a") (logy), 

and, therefore, for c = 1 — 1/(2 logy), that 

ip{z,L 2 )< 6.4193 ■ z a • .3764 (log x) 2/3 log log x 

< 2.4163 ■ x° ■ T (log log x) 1+a ■ (logx) 2/3 . 

Thus, we want to find C\ such that 

4.8325 • x ■ (log log a;) 2 • (loga;) 2/3 < x 

^l/(.7528(log a: ) 2 / 3 loglogx) - 1 I' 

Plugging in xo, we may take C\ — 13553617.97. 

It turns out that the other constants are negligible compared to C\ , so we will 
sacrifice some sharpness in bounds for space in what follows. 



Ci'- Note that we can bound 



Urn / j \ 1,2 1,3 } Z. ^ U2 jU 

fcm>/3 K 3 \ K K / 3 K K 

k - k>l2 k>t2 

m>2 

V- 1 1 nnm 1 1-OOOlz 

= -E^T)<i-oooi,E^<7kT 

for x > xo- Now we need 

1.00012 „ x 
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Plugging in x = xq into the formula for I shows that we may choose Ci 
.1862. 



C3 : A simple application of Lemma 14.11 gives 

x 1 x ( 4 + 3 log log (x/p) 

P .Q ~ P\ P 

g<a:/p 

Now, summing this over p > d A gives at most 

,» o, , \ 1 ^ (4 + 3 log log x) 
x(4 + 3 log log x)J2 "a < ■ 



by Therefore, we need 



.t (4 + 3 log log x) C3X 
£ 4 - 1 " ~T~' 

which holds for CS > 4.3 x 1CT 12 . 



C4: Here, we need 

4a; log log x ^ C42; 



L 2 £ 
4 log log 2; exp ^(logx) 1 ^ 6 



< C* 4 . 



exp ^.61 (logx) 2 ^ 3 log log 2;^ 
Plugging in x = x gives that we may choose any C4 > 10~ 17471 . 

C5: We must next compute C5 such that the inequality 

4e ln225 x log logic (1 + 21oglogx) ^ 1 < C$x 



E 



p ^ P - t 

P <x/e 1 



holds for all x > xq. By ([5]), this requires 
4 e 1.0225 j og i g x (1 + 2 log log x) 
exp ((loga;) 1/6 ) 

For x > xq, we may take C5 = 2.8117. 



(log log ie + .2615) < C 5 . 
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Cq: We have 

E E E E ^^-EEEsE?- (17) 

r>£ 4 <j=-l(r) m=0(g) g'=-l(r) r q m q> 

q<x m<x <?'<z 

Applying Lemma I47T1 we find (|17|) is bounded by 

z(4 + 31oglogz)^^^l. (18) 

r q m 



Now, since \og(x/q) < log a;, we can bound 

E 1= E - = 1 E 

m=0(c, 



1 1 v-^ 1 logic 
m nq q n q 

\){q) n<x / q n<x/q 



Using this and another application of Lemma 14.11 the iterated sum in (fl8|) 
above can be bounded by 

z(4 + 3 log log 2;) (log x) (4 + 3 log log x)^ \. (19) 

r 

Lastly, (|12l) shows that 



V — 1 



and so, putting this into (|19p . we need 

z (4 + 3 log log z) (log sc) (4 + 3 log log x) x 

—1 -° 6 T (20) 

A calculation shows that (J20J) holds for C 6 = .0054. 

Thus, for x > exp(l0 6 ), the number of amicable numbers not greater than x, 
A(x), can be bounded as 

A{x) < 13553620.97 



exp ^(logx) 1 ^ 6 



Then, by partial summation, we have that the third sum in ([3]) is bounded by 

13553620.97 (exp (-10) + 1) / -. r- < 654666169. 

Jxo texpf(logi) 1/6 J 

Combined with the results of sections [2] and [3l this proves Theorem 11.11 
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